Generalized entropy production fluctuation theorems 

for quantum systems 

Shubhashis Rana* Sourabh Lahiri^ and A. M. Jayannavar-f 

Institute of Physics, Sachivalaya Marg, Bhubaneswar - 751005, India 

Abstract 

Based on trajectory dependent path probability formalism in state space, we 
derive generalized entropy production fluctuation relations for a quantum system 
in the presence of measurement and feedback. We have obtained these results for 
three different cases: (i) the system is evolving in isolation from its surroundings; 
(ii) the system being weakly coupled to a heat bath; and (iii) system in contact with 
reservoir using quantum Crooks fluctuation theorem. In case (iii), we build on the 
treatment carried out in [H. T. Quan and H. Dong, arxiv/cond-mat: 0812.4955], 
where a quantum trajectory has been defined as a sequence of alternating work and 
heat steps. The obtained entropy production fluctuation theorems retain the same 
form as in the classical case. The inequality of second law of thermodynamics gets 
modified in the presence of information. These fluctuation theorems are robust 
against intermediate measurements of any observable performed with respect to 
von Neumann projective measurements as well as weak or positive operator valued 
measurements. 

PACS: 05.40.Ca, 05.70.Ln, 03.65.Ta 

1 Introduction 

Nonequlibrium processes are common in nature, but a general framework to understand 
them is lacking as compared to equilibrium systems. However, recent development in 
the field of nonequilibrium statistical mechanics, had led for the discovery of fluctuation 
theorems (FT) [Ij-tZl, which are exact equalities that are valid even when the system of 
interest is driven far away from equilibrium. For such a nonequlibrium system, the sta- 
tistical distribution of thermodynamic quantities such as work, heat, entropy etc. exhibit 
universal relations. These thermodynamic quantities have now been generalized to a sin- 
gle trajectory of system evolving in phase space. They are random variables depending 
on the phase space trajectory (stochastic thermodynamics [S]). The physical origin of 
FTs rely on the time reversal symmetry of the dynamics and they are expected to 
have important applications in nanoscience and biophysics. The second law of thermo- 
dynamics emerges in the form of inequalities from these theorems [3l[6] . It can be shown 
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that second law is valid on average. Here averaging is done over different trajectories, 
thus not ruling out the possibility of transient violations of second law for individual 
realisation [9]. These theorems have helped us in understanding how thermodynamic 
irreversibility arises from the underlying time reversible dynamics |1U] . 

One of the FT was initially put forward by Jarzynski |3] in the form of the nonequi- 
librium work theorem, by means of which one can extract information about equilibrium 
changes in free energy AF by measuring the nonequilibrium work W performed on a 
system by the external drive. The system is initially prepared in equilibrium, and then 
driven away from equilibrium using some predetermined protocol X{t) which runs from 
t = to t = T. The Jarzynski Equality is given by 

(e-^W) = e-^^^. (1.1) 

The work W, depends on trajectory of the system, whose initial state is sampled from 
equilibrium distribution. The angular brackets denote averaging over an ensemble of 
such trajectories and the free energy differences AF = F{X{t)) — F{X{0). A stronger 
fluctuation theorem was provided by Crooks in the form 



Pri-W) 



^Pi.w-^,)^ (1.2) 



Pf(W) and Pr{W) being the work probability densities generated under the forward 
protocol \(t) and the reverse protocol A(r — t), respectively. 

A more general FT was put forward by Seifert [7] which contains the Jarzynski and 
the Crooks theorems as special cases. A system which is in contact with a heat bath, 
is initially prepared in some arbitrary distribution ^0(3^0) of phase space points and is 
perturbed by varying an external parameter X{t) up to time t = t. In the reverse process, 
the system evolves from some other initial distribution ^1(0;,-) under the time-reversed 
protocol A(r — t). The Seifert's fluctuation theorem states that, the probability of a phase 
space trajectory along the forward process, P[x(t)] is related to that along the reverse 
process, P[x{t)], as 

P[x{t)] P[xit)\Xr]pi{Xr) VlM 

where ASb is the change in the entropy of the bath {ASb = ^, where Q is heat absorbed 
by the bath). P[x(t)|xo] is a short notation for functional of a path starting at xq and x{t) 
being the phase space trajectory ending at Xr- If, in particular, the distribution pi{xr) 
in the final distribution at time r as dictated by dynamics, then the above relation gives 
the integral fluctuation theorem (IFT) for total entropy production [7|: 

(e-^^'°*) = 1. (1.4) 

where 

AStot = AS + ASb = In ^ + ^. (1.5) 



Po(^'o) 

refer to Seifert's article |7]. If the system is in steady state, one can also obtain detailed 



Here In is the change of system entropy along a given trajectory. For details we 
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entropy production fluctuation theorem (DFT), namely 

e^^"". (1.6) 



The IFT follows directly from the DFT. Using Jensen's inequality in Eq.f ll.4p we get 

(AStot) > 0. (1.7) 

This is a statement of second law of thermodynamics, expressed in the form of inequality 
for the average change in total entropy. 

If the systems are driven by the feedback controlled protocols, which in turn depend on 
the measurement outcomes of the state of the system at intermediate times (information 
gain), then IFT gets modified a form [TT] 

^^-AStot-I^ = 1, (1.8) 

where I is the mutual information which quantifies the change in uncertainty of state of 
the system upon making measurements. Application of Jensen's inequality generalizes 
second law for total entropy production: 

{AStot) > -(/). (1.9) 

The average mutual information (J), is always non negative [12]. Thus the average 
entropy change can be made negative by feedback control, and the lower bound is given 
by — (/). There are few attempts to extend the IFT (Eq. fll.4p ) to the quantum domain 
[TSHIS]. In our present work we extend the IFT for AStot to quantum systems, in presence 
of multiple measurements and feedback. We assume that measurement procedure involves 
errors that are classical in nature. We show the robustness of FTs against intermediate 
measurements of any system observable (both von Neumann projective measurements or 
generalized positive operator valued measurements (POVM)). 

We obtain these theorems for three different cases: (i) the system evolves in isolation 
from its surroundings; (ii) it is weakly coupled to a heat bath; and (iii) evolution of system 
coupled to heat bath is modelled in terms of work steps and heat steps following closely 
the treatment given Ref (TB] and is described in section HJ Our treatment is based on 
path probability in state space. The measurement is assumed to be von Neumann type, 
i.e, projective measurement which results in the collapse of system state to one of the 
eigenstates of the corresponding observable. Case (i) namely isolated quantum system 
is discussed in detail. DFT is obtained for various different situations, i.e., (a) system 
evolving unitarily, (b) in the presence of measurement and feedback, and finally (c) in the 
presence of intermediate measurements, of any observables of the system. The IFT follows 
from DFT. For cases (ii) and (iii) we have derived generalized IFT. In the appendix, we 
have given a proof of IFT in presence of weak measurements. In passing, we note that 
all the extended quantum FTs retain the same form as their classical counterparts. 



2 Isolated quantum system 
2.1 Unitary evolution 

In this section we consider an isolated quantum system given by Hamiltonian H{X{t)), 
where A(t) is some external time dependent protocol. To clarify our notation and for 
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completeness we rederive DFT for this system following the treatment of ref [H]. Initially 
at time t=0 energy measurement is performed and system is found to be in eigenstate 
\io), with energy eigenvalue Eo. It then evolves unitarily from time to r under the 
protocol A(t). The energy measurement at final time r is performed and system is found 
to be at state |v) with energy eigenvalue E^-. If the initial probability density of the state 
l^o) is p{io) then the joint probability of l^o) and |v) (forward state trajectory) is given 
by 

PF{iT,io) = p(vKo)p(^o) 

= |(v|f/A(r,0)Ko)lV^o), (2.1) 

where U\{t2, ti) denotes the unitary evolution operator for given A(t) from time ti to time 
t2- It is defined as 

Uxih, tl) = Texp (^-^ ^ ' H{X{t))d?j . (2.2) 

Here, T denotes time ordering. 

The system entropy is defined as S(t) = — In p{it). As the system is isolated there is 
no generation of heat, i.e, Q=0, Using eq. (11. 5p the total change in entropy production 
AStot during the evolution from time to r is equal to the change in system entropy 
alone. 

AStot = -\nP^, (2.3) 

where p{ir) is the final probability of state |v) at time r. The probability density 
PpiAStot) for the forward path is by definition 

PpiAStot) = ^ (^^tot + In Pf(v, ^o) 

= {"^^^ot + 1^^) Pi^r\^o)p{^o). (2.4) 



We now introduce time reversal operator G. The time reversed state of \i) is defined as 
\i) = Q\i). It can be readily shown that [TT] 

p{i2\ti) = \{t2\Ux{t2,ti)\ti)\^ = {{tilU^tihMM^ = PiH\~^2). (2.5) 

where t = t — t and A^(t) = A(r — t) is the time reversed protocol of A(t). The evolution 
of the system from given time reversed state Q\i2) to the time-reversed state 0|ii), under 
the time reversed protocol A^(t), is given by the conditional probability p{ii\i2)- We 
consider the initial distribution of reverse trajectory to be equal to the final distribution 
of forward trajectory 

P{l)=p{tr). (2.6) 

The states \i) and \i) have one-to-one correspondence. Multiplying and dividing by pi^ir) 
in the summand in eq. fl2.4p and using (12. 5p and (12.60 . we get 

PpiAStot) = Y^S ( AStot + \n^] p{^o\^r)p{^r)^ 

rf V pW)/ p(v) 



It, 10 
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7- in ^ PVt) J 



AStot 



Pni-AStot)- (2.7) 

To arrive at this result we have used eg. (12. 3 p in the second step and eg. (12.61) in third step. 
-Pr(v,^o) is the joint probabihty of the corresponding states in the reverse direction. If 
the AStot is the total entropy change for forward path then the total entropy change in 
the corresponding reverse path is —AStot- It follows from the fact that p(v) and p(^o) is 
the initial and final probability distribution of state in the time reversed process because 
of unitary evolution. The eg. (l2.7p can be written in the form 

- e^^-. (2.8) 



Pni-AStot) 



This is the detailed fluctuation theorem for change in total entropy, extended to the 
guantum regime. Simple cross multiplication followed by integration over AStot leads to 
the integral form of the above theorem: 



g-A5t.t^ = 1. (2.9) 



2.2 Isolated quantum system with feedback 

So far we have been dealing with a predetermined protocol, also known as open loop 
feedback. Often to increase the efficiency of a physical process ( eg. engines at nanoscale, 
molecular motor etc.), we need to perform intermediate measurements and change the 
protocol as per the outcomes of these measurements pTH2l] . Such a process is known 
as closed loop feedback. Let the system evolve under some external protocol Ao(t), from 
its initial energy eigenstate \io) measured at time to- At time ti, we perform a measure- 
ment of some arbitrary observable and system collapses to state We assume that 
measurement process leading to information gain involves classical errors. Here yi is the 
measured outcome with a probability p{yi\ii), while system's actual state is De- 
pending on the value of yi the protocol is changed to \yi(t). Under this new protocol the 
system evolves unitarily up to time t2 where another measurement is performed and so 
on. This process terminates at time r when the system collapses to its final energy eigen- 
value We should note that initial and final measurements are energy measurements. 
The joint probability of the corresponding state trajectory for n number of intermediate 
measurements yi,y2, ■ ■ ■ Vn at times ti,t2,- ■ - tn respectively is [11] 

Pf(v, ..,ii,io,yn, ..,yi) = p(vKn) ■ ■ ■ piy2\i2)pii2\ii)piyi\ii)piii\io)p{io) 

= |(v|t/A,„(r,ti)|Op...p(y2K2)|(^2|f/A,,(t2,tl)Kl)lMz/lKl) 

X \{^l\Ux,{h,0)\^o)\M^o)■ (2.10) 
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It may be noted that the joint probabihty of path is expressed using classical probability 
rules. This is because we perform projective measurement on the system which col- 
lapses to one of the eigenvalues of the measured observables [TT1|25]. As a consequence, 
it wipes out the previous memory of evolution and the post-measurement evolution be- 
comes uncorrelated to the pre-measurement evolution. Thus if one performs intermediate 
measurements along two paths, the interference effects between the two paths disappear 
and the quantum effects are suppressed. Hence in the presence of measurement, path 
probability in state space obeys classical probability rules, and is given by product of tran- 
sition probability of paths between consecutive measurements. However, it may be noted 
that quantum mechanics enters through the explicit calculation of transition probability 
between two consecutive states. 

To generate the reverse trajectory of a path in state space given in eq.f l2.10p . we first 
choose one of the forward protocols with probability p{yn ■ ■ ■ ,1/2,1/1), and then blindly 
time reverse the protocol. We perform measurements at the appropriate times along 
reverse path to allow the state to collapse to the corresponding time-reversed eigenstates. 
We do not use these measurements to perform any feedback to respect causality [23] . 
Then the expression for the joint probability of reverse trajectory is given by 

Pniir--- ,io,yn,--,yi) = p(«7v|v) ■ ■ ■p(^oKi)p(v)p(l/n ■ ■ ■ ,yi)- (2.11) 

The mutual information gain due to measurements between the measured values and the 
actual value is defined as pT|[23] 

/ = i^ P^yM--Piy2\i2)p{yi\ii) .2 12) 

p{yn--- ,y2,yi) 

We now calculate the joint probability density Pp^AStot,^) of the entropy production 
and X along the forward path, which is 

PpiAStot,^) = j dun- ■ -dyi ^ 5 ^ASfot + In 5 (X - /(i^, .., zi, ?/„, .., 2/1)) 

X PF{ir,--,il,io,yn,--,yi) 

dyi ^ (5 fAStot + ln^^|^ j (5(X--/(i„,..,2i,?/„, ..,1/1)) 

ir---,h,io ^ y\ J / 

X p{ir\in) ■ ■ ■p{y2\i2)p{i2\ii)p{yi\ii)p{ii\io)p{io) 

/ p(i )\ 

dyi 22 ^ i^Stot + \n—^j 6{I - I{in,..,ii,yn,--,yi)) 

iT---,h,io I / 

X p(iAr|v) ■ ■ •p(«oKl)p(v)p(l/n " " " , l/l) 6^^*°*+^ 

dyx 2I1 (A'S'tot + In— 5 (X- /(«„,. .,«!,?/„, ..,1/1)) 

X Pr(%^--- ,Zo,2/n, -,2/1)6^^*°'+^ 

jdyn---dyi ^ 5 |^AS'tot + In 5 (X - /(i„, .., ii, ?/„, .., ?/i)) 

X Pniir ■ ■ ■ ,io,yn, ■■,yi) 



dy„ 



dy„ 



dy„ 
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= e^^-'+^P^(-A5tot,X). (2.13) 

In deriving above result we have made use eq. fl2.10p . (12. lip . (I2.12p . The path variable 
J(z„, .., ii, .., ?/i) is given by eq. (I2.12p . and X denotes its value. It is important to 
note that the probability density function Pr(— AS'tot,X) gives the probability of reverse 
trajectories along which the entropy chage is —AStot and whose corresponding forward 
trajectory has the mutual information X between its measured outcomes and actual states. 
Once again, the initial and final distributions of states along forward trajectory get in- 
terchanged in the reverse trajectory because of unitary evolution between measurements. 
Along the reverse trajectory the change in total entropy is —AStot- Thus we obtain the 
DFT 

PpjAStot,^) _ ^AStot+x_ (^2.14) 



Pr{- AStot,!) 

From the above equation the extended version of IFT and second law, eqs. (II. 8p and 
(II. 9p can be readily obtained as discussed in earlier subsection. 

2.3 Isolated system under multiple measurements 

In this subsection we restrict ourselves on the influence of intermediate measurements 
of arbitrary observables on the statistics of AStot- To this end we do not involve any 
feedback. Following closely the discussions in section (12. 2p . of path probability in state 
space is given by 

P{ir, --,ii,io) = piirlin) ■ ■ ■ p{i2\ii)piii\io)p{io) 

= |(v|?7A,Jr,ti)|z„)|2...|(^2|f/,^^(t2,ti)Ki)n(^i|f/Ao(ti,0)|zo)rp(^o). (2.15) 

From preceding section we now calculate the probability density PpiAStot) of the total 
entropy change along forward path 



"Tl --1 "il 



P,{AStot)= E 6(AStot + \n^)p{ 

= 2J 5 ( AStot + In— ^ j p(v|z„) ■ ■ ■p(z2|ii)p(ii|zo)p(^o) 

It— ,il,i0 y\ J / 

= 5(A5tot + ln^')p(z;v|v)---p(^o|n)p(v)e^^- 

J2 ^ {^^Stot + In ^) Pnil---, ^o)e^^'- 
- V S ( AStot + \n^] Pnil- ■■ ;zo), (2.16) 



It--- ,«1,«0 

It--- ,il,io 



where Pniir ■ ■ ■ , ^o) is the probability of reverse path. The DFT for AStot follows from 
the above equation: 

PpjAStot) _ AStot (n 1 7\ 

PRhAStot) ~ • ^ ^ 
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We observe from eq. fl2.17p the robustness of this FT against intermediate measurements 
[26l|27]. It retains the same form as in the classical case. The path probability, however, 
gets modified in presence of measurements and statistics of AStot is strongly influenced by 
the intermediate measurements. In the next section we derive IFT in presence of feedback 
for a quantum system coupled weakly to a bath. In the appendix, we have shown that 
the IFT for AStot is also robust against weak or generalized intermediate measurements. 



3 Weakly coupled quantum system 

Consider a driven system which is weakly coupled to a bath. The total Hamiltonian will 
be 

HiXit)) = Hsim) + Hb + Hsb. (3.1) 

The external time dependent drive X{t) only affects the system Hamiltonian Hs{X{t)) 
while the bath Hamiltonian Hb and interaction Hamiltonian Hsb are time independent. 
As the system is weakly coupled it is assumed that Hsb is negligibly small compared to 
Hs{X(t)) and Hb- Initially the super-system (system+bath) is coupled to a large reservoir 
of inverse temperature P |27l|28] . At time t = the large reservoir is decoupled from the 
super-system. Hence initially the super-system will be in a canonical distribution, 

-/3/f(A(0)) 

where y(A(0)) = Tr e"'''^'-^*^'^^-' . The system and the bath Hamiltonians commute with 

each other, hence we can measure simultaneously the energy eigenstates for system as 
well as bath. At t=0, the measured energy eigenvalues of system and bath are denoted by 
Eq and Eq , respectively. We perform N number of intermediate measurements of some 
arbitrary observable at time ti,t2- ■ - tN between time to r. Initially the protocol was 
Ao(t). At ti the measured output is yi, while its actual state is ii, with probability 
p{yi\ii). Now the protocol is changed to Aj,^(t) and system evolves up to time t2 ■ 
Again measurement is performed and protocol is changed according to the output at 
intermediate times and so on. Finally at t=r joint measurement is performed on system 
and bath Hamiltonians, and the measured eigenvalues are E^ and E!^, respectively. The 
system-reservoir interaction energy can be neglected in the presence of weak coupling. 
Hence during the evolution process from time t = to t = t for a single realization the 
change in the internal energy of the system is given by [28] 

AU = E^ - E^ (3.3) 

and the heat dissipated to the bath is 

Q = E^- <. (3.4) 

If io and v denote initial and final system energy eigenstates, then system entropy change 
is 

ASsys = -ln^, (3.5) 
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and the total entropy change is 



AStot = AS + = - In ^ + |, (3.6) 

where T is the temperature of the bath. The mutual information between the state 
trajectory {^l, ^2, ■ ■ ■ , "^tv} and the measurement trajectory {yi, ?/2, ■ ' ' Un} is 



J = ln 



pivilh) ■ ■ -piyNliN) 



(3.7) 



Denoting initial and final states of the bath by and a,-, it can written from microscopic 
reversibility pT|l2i] 

p{ir, OrNo, tto) = P{io, "olv, "r)- (3.8) 

where p{iT, OtI^O; ^o) is the total transition probability for system and reservoir to evolve 
from state \io,ao) to \ir,(yr) under the full Hamiltonian. Here \i,a) = Q\i,a) is the 
time-reversed state of \i,a). To generate the reverse trajectory, proper causal protocol 
has to be used which has been discussed in section 12.21 Thus the forward and the reverse 
path probabilities of trajectories are respectively given by 

Pf{A ^ B)= p{ir, ttrliN, On) " " " p(2/i Ki)p(^i, "iKo, o:o)p{io, ao), (3.9) 

Pr{A ^ B)= p{io, do\ii, ai) ■ ■ -piiN, "ivlv, ar)p(v, ar)p(l/i, 2/2, ■■ ■ Vn)- (3.10) 

The notations A and B denote initial and final values of the forward protocol, respec- 
tively. For reverse trajectory we have chosen the outcome of the forward trajectory with 
probability p{yi,y2, ■ ■ - Vn) and have blindly reversed the protocol, but performing mea- 
surements (without any feedback) at appropriate time instants. From fl3.9p and fl3.10p 
we get 

Pf{A -^B) _ 
Pn{A ^B)~ 

^ p{yN\iN) ■ ■ -piyilH) pjio^ao) 
PiVi, ■■■ ,yN) piir,ar) 
^ p(^o)p(Qo) 

p{ir)p{ar) 

(3.11) 

In arriving at (13. lip , we have used microreversibility (13. 8p and we have assumed that the 
system and the bath are weakly coupled. The joint probability of system and bath states is 
approximated as a product of individual state probabihties. Correction to this factorized 
initial state is at least of second order in system-bath interaction, and therefore they 
can be neglected in the limit of weak coupling. The bath probability can be considered 
canonical with inverse temperature /3. This leads to 

Pf{A B) _ J AS.„, /Zb _ I ASsyspQ/T _ AStot+I /o ion 

Pn{A ^B)-' ' e-^^?/Z^ - e e e - e . (3.12) 



p{ir,Olr 


iN^oiN) ■ ■ -piyi 


|ii)p(ii,ai|2o,ao)p(^o,ao) 


p{io,ao\ii,ai) ■ ■ ■p{iN,aN\ 


ir,ar)p{ir, 


ttr)p(2/i,2/2, ■■■Vn) 
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A simple cross multiplication and integration over paths gives the extended IFT. It may 
be noted that in our framework we can also obtain the DFT, provided the system either 
begins and ends in equilibrium or in the same nonequilibrium steady state [5]. In the 
next section, we prove the same IFT for AStot by means of the method developed in [16j 
by using the quantum mechanical generalization of the Crooks fluctuation theorem. 

4 IFT using quantum Crooks fluctuation theorem 

We consider the system to be coupled to a bath, but there is no assumption made in 
regard to the strength of the coupling. Each time step in the entire evolution is divided 
into two substeps. In first substep the protocol is changed while in second, protocol is 
kept fixed and system relaxes by dissipation of heat. The total evolution is divided into 
N steps. Each step starts at t„ and ends at t^+i, where n = 0, 1, 2, ■ ■ ■ — 1. We closely 
follow the treatment in [16] . 

For a quantum adiabatic process the protocol changes slowly and the system remains 
in same eigenstate in the work step. However, in the present case the work step is almost 
instantaneous and the process is non adiabatic. As a consequence the eigenstates before 
and after work step may be different. The system starts to evolve under a predetermined 
protocol Aq. For simplicity, let us consider the observable measured at intermediate times 
to be the Hamiltonian itself. It can be readily generalized to the case of other observables. 
We consider that the feedback is applied at the beginning of the each work step and we 
change the protocol subsequently according to the result obtained from the measurement, 
as discussed earlier. The conditional probability p{yn-i\in~i) denotes that the measured 
outcome is i/n-i while the actual collapsed state is \in-i, An-i), at the beginning of the n*'^ 
work step. Within the ket notation, i„_i represents the state of the system and A„_i is the 
value of the control parameter. After the measurement of tn-i, the protocol is changed to 
An(?/n-i) from A„_i(?/„_2). During the work step, the system evolves unitarily from 
to t'^_i, where it is measured to be in state A„). The time taken in the work substep 
is considered to be too small for the system to relax. In the n*^ heat step, the system 
relaxes from state A„) to \in, A„). Therefore, the path followed by the system in state 
space of the measured eigenstates from state \io, Xq = A) to Iv^A,- = B) is represented 

as \io,Xo) 1^0, Ai) |i'i,A2) > \iN-i,XN~i) -> \i'N-i,^N) \iN,XN)- 

Let E{in, An) be the energy eigenvalue of state \in, A„). By adding the contributions from 
all the work steps, the total work done on the system is given by 



N-l 




(4.1) 



n=0 



while heat dissipated into the bath is 




(4.2) 



n=0 



The change in internal energy of the system along the trajectory is 



AE = Q + W = E{iM, \n) - E{io, Ao). 



(4.3) 
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As before, the mutual information is 

J ^ p{yn\in)-p{y2\i2)p{yi\ii) ,^ ^^ 

Pivn--- ,y2,yi) 

The forward and the reverse path probabihties are respectively given by 

7V-1 

Pf{A ^ B)= p{io, Ao) Yl Piyn\in)PF{\in, A„) \n+l)) PfHC, A„+i) \in+l, Ki+l)) ■ 

n=0 

and 



(4.5) 



(4.6) 



n=0 



As mentioned earlier, during the work step, the system can be regarded as an isolated 
quantum system and evolution is completely determined by the time-dependent Hamil- 
tonian Hs{X{t)). Thus the evolution is unitary. Microscopic reversibility for work step 
gives [16] 

Prilin, A„) Xn+i)) = PRiK, A„) ^ A„+i)). (4.7) 

The heat steps or relaxation steps are assumed to be microscopically reversible and 
obey the local detailed balance for all the fixed values of the external parameter A. The 
detailed balance condition in relaxation substep implies 

-Pf(K) A„+i) -> |i„+i, A„+i)) r X 77/-' \ M fA Q\ 

p X ^ T ^ ^ = exp -/3(K+i, A„+i - E{t^, Xn+i)]- (4.8) 

Using above two equations we get 

Pf{A B) p{yn\in)--p{y2\i2)p{yi\ii) p(^o,Ao) j^c \ \\ 

p f A ^ m = 7 ^ r — 11 exp -pO^n+i, A„+i - E{i^, Xn+i) ■ 

Pr{A ^ B) p[yn ■ • ■ , 1/2, yi) p{^N, A^v) Y 

(4.9) 

The total entropy change along the trajectory, AStot = AS + ASb, which is a trajectory 
dependent random variable, where AS" = — In '^^(^'j^^^'* is the change in system entropy, 
and ASb = Q/T is the entropy change of the bath, along a single trajectory. Using (14. 2 p 
and f Ojl . eq. (K9\i simplifies to 

PpjA B) _ J AS,y Q/T _ AStot+I (A 1 n\ 

P^{A^B) ^ ' 

This immediately leads to the generalized integral fluctuation theorem for total entropy 
change, in the presence of measurement and feedback. In the above derivation, we have 
taken measurements for feedback at the beginning of the work steps for simplicity. These 
measurements can be performed at any time in-between the work steps. The result will 
not be affected. It would only make the notations more complicated and would not 
provide any new physical insight. Feedback cannot be performed within the heat step 
which be definition requires protocol to be held constant. 

As in case (ii), the DFT for AStot can be obtained if the initial and final distributions 
are in the equilibrium or in the same nonequilibrium steady state |5j. 
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5 Conclusions 



Based on the path probabihty formulation in state space, we have derived generahzed total 
entropy production fluctuation theorems for quantum systems in presence of measurement 
and feedback, for three different cases. They retain the same form as in classical case. 
The second law of thermodynamics gets modified in the presence of information and 
feedback (eq. (11.91) ). For isolated quantum system with feedback, we have derived the 
generalized DFT for the total entropy. For this case DFT retains the same form in 
presence of multiple measurements of any system observable, thus showing the robustness 
of these fluctuation theorems against measurements (von Neumann type or generalized 
measurements). For the case (ii) of a weakly coupled quantum system under feedback, 
we have derived the extended IFT for total entropy. In case (iii), we have derived the 
extended IFT for AStot, using the quantum Crooks fluctuation theorem, where quantum 
trajectory is characterized by a sequence of alternating work and heat steps. IFT is 
valid for any initial arbitrary state of a system. DFT in cases (ii) and (iii) can be 
obtained only when the system either begins or ends in equilibrium or remains in the 
same nonequilibrium steady state. By using our approach, the generalized DFT can 
be proved, but we have not provided the details. The derivation of the robustness of 
the fluctuation theorems against intermediate measurements is given only for case (i), 
namely, for the isolated quantum system. Following the same treatment, the robustness 
of fluctuation theorems can be readily demonstrated for cases (ii) and (iii) as well. 

In conclusion, we have generalized total entropy production fluctuation theorem in 
presence of feedback to the quantum domain using three different approaches. 
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Appendix 



A Isolated system under weak measurements 

In this appendix we derive IFT for AStot under weak measurement (POVM), as opposed 
to projective von Neumann type measurements considered in sec ]2.1[ We follow the math- 
ametical treatment given in [29j. For simplicity we consider only one weak measurement 
is performed at intermediate time. The genaralization to multiple weak measurements is 
straight forward. Consider an isolated quantum system is controlled externally through 
time dependent protocol A(t). Initially at t = energy measurement is performed and 
the system is found to be in state \n, 0) with probability density p„. The density matrix 
becomes 

P„(0+) = ^^^^ (A.l) 

Pn 

where po is the density matrix of the system before measurement and 11° denotes von 
Neumann projective measurement operator and p„ = Tr n°po- The system then evolves 
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unitarily up to time ti and a weak measurement is performed and we get the density 
matrix 

Pnitt) = 5^M,?7A(ti,0)p„(0+)t/l(ti,0)Mj. (A.2) 

r 

Mr is the weak measurement operator with property Y^^M^M^ = 1. The system un- 
dergoes further unitarily evolution and finally the projective measurement is performed 
and the system is found to be in state |m, r) at time r. The conditional probability 
P\{m, T\n, 0) for system initially in state \n, 0) and finally in state \m, t) is given by 

PAKr|n,0) = Tr IiljJ^{T,h)p^{t+)Ui{T,h). (A.3) 
Thus the probability of the change in total entropy, 

^Stot = - hipm + Inpn, (A.4) 
for a given pre-determined protocol \{t) is 

Px{AStot) = ^S{AStot + ^^Pm-'^^Pn)P\im,T\n,0)pn. (A.5) 

m,n 

where pm is the probability of the system to stay at the end of protocol at final time r. 
The Fourier Transform of this probability is 

Gx{u) = j d/\StotPx{AStoM^''''''\ (A.6) 



Substituting the expression for Px{AStot) from eg. flA.51) and using eqs. flA.3p . flA.2p . flA.ip 

we get 

m,n,r 

= ^ Tr n^e-»'"^/f/A(r,ti)M,?7A(ti,0)nOe^"^'^^Vof/l(ti,0)Mtf/t(r,ti) 

m,n,r 

= J]Tr e-^'^i°''/f/A(r,ti)M,f/A(ti,0)e^"^^''Vof/l(ti,0)Mt?7l(r,ti) 



5^Tr [/J(ti,0)Mtf/t(r,ti)e-^"'"''^t/,(r,ti)M,f/A(ti,0)e»i'^^Vo 



(A.7) 



p/ is the final density matrix which is diagonal in the energy basis. In the second step 
we have used the completeness relation J2m = 1 J2n = ^ ^^e projective 
operator. Substituting u = i, for the Fourier Transform variable, we get from eg. (1A.6|) 

Gx{i) = (e"^'^'°*) and hence 



r 

= J]Tr Ul{tuO)M}ul{T,h)e'^PfUx{T,U)MrUx{tuO) 



Inpo 



Po 
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= J]Tr Ux{r,t,)MrUx{h,0)Ul{h,0)Mlul{T,t,)e'^Pf 

r 

= ^Tr UxiT,t,)Mr.Mlui{T,t^)e'-''f 

r 

= Ti:Ux{T,h)Ui{T,h)e'^Pf 

= Tr 6'""^ =Tt pf = 1. (A.8) 

In the second line , we make use of e^^'^^^po = 1? while the cyclic property of trace is 
used in the third line, operator identity MrM^ = 1 is used in fourth step. 

We have proved that the IFT holds in the same form as in eq. (11.41) under intermediate 
weak measurements. This can be readily generalized to the multiple intermediate weak 
measurements, which corroborates the robustness of fluctuation theorems under weak or 
generalized measurements. 

References 

[1] D. J. Evans, E. G. D. Cohen and G. P. Morriss, Phys. Rev. Lett. 71, 2401 (1993). 

[2] D. J. Evans and D. J. Searles, Phys. Rev. E 50, 1645 (1994). 

[3] C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997). 

[4] G. E. Crooks, J. Stat. Phys. 90, 1481 (1998). 

[5] G. E. Crooks, Phys. Rev. E 60, 2721 (1999). 

[6] U. Seifert, Phys. Rev. Lett. 95, 040602 (2005). 

[7] U. Seifert, Eur. Phys. J. B 64, 423 (2008). 

[8] K. Sekimoto, Prog. Theor. Phys. Suppl. 130, 17 (1998). 

[9] M. Sahoo, S. Lahiri and A. M. Jayannavar, J. Phys. A: Math. Theor. 44, 205001 
(2011). 

[10] R. Kawai, J. M. R. Parrondo, and C. Van den Broeck, Phys. Rev. Lett. 
98, 080602 (2007). 

[11] S. Lahiri, S. Rana and A. M. Jayannavar, J. Phys. A 45, 065002 (2012). 

[12] T. M. Cover and J. A. Thomas (2006) Elements of Information Theory, 2nd edn 
(Hoboken, NJ: Wiley-Interscience) 

[13] Massimiliano Esposito and Shaul Mukamel, Phys. Rev. E 73, 046129 (2006). 

[14] T. Monnai, Phys. Rev. E 72, 027102 (2005). 

[15] S. Deffner and E. Lutz, arXiv/cond-mat: 1201.3888. 

[16] H. T Quan and H. Dong, arxiv/cond-mat: 0812.4955. 



14 



[17] A. Hopkins, K. Jacobs, S. Habib and K. Schwab, Phys. Rev. B 68, 235328 (2003). 

[18] D. Steck, K. Jacobs, H. Mabuchi, T. Bhattacharya and S. Habib, Phys. Rev. Lett. 
92, 223004 (2004). 



[19] K. Jacobs, arxiv: quant-ph/0605015 



[20] K. Jacobs and A. Shabani, Contemp. Phys. 49, 435 (2008). 

[21] T. Sagawa and M. Ueda, Phys. Rev. Lett. 104, 090602 (2010). 

[22] T. Sagawa and M. Ueda, Phys. Rev. E 85, 021104 (2012) 

[23] M. Ponmurugan, Phys. Rev. E 82, 031129 (2010). 

[24] J. M. Horowitz and S. Vaikuntanthan, Phys. Rev. E 82, 061120 (2010). 

[25] S. Rana, S. Lahiri and A. M. Jayannavar, arxiv/cond-mat: 1111.6229. 

[26] M. Campisi, P. Talkner and P. Hanggi, Phys. Rev. Lett. 105, 140601 (2010). 

[27] M. Campisi, P. Hanggi and P. Talkner, Rev. Mod. Phys. 83, 771 (2011). 

[28] P. Talkner, M. Campisi and P. Hanggi, J. Stat. Mech. , P02025 (2009). 

[29] M. Campisi, P. Talkner and P. Hanggi, Phys. Rev. E 83, 041114 (2011). 



15 



